Electronic Structure of Sodium Cobalt Oxide: Comparing Mono- and Bilayer-hydrate 
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To shed new light on the mechanism of superconductivity in sodium cobalt oxide bilayer-hydrate 
(BLH), we perform a density functional calculation with full structure optimization for BLH and 
its related nonsuperconducting phase, monolayer hydrate (MLH). We find that these hydrates have 
similar band structures, but a notable difference can be seen in the a lg band around the Fermi level. 
While its dispersion in the z direction is negligibly small for BLH, it is of the order of 0.1 eV for 
MLH. This result implies that the three dimensional feature of the ai g band may be the origin for 
the absence of superconductivity in MLH. 
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Layered sodium cobalt oxide bilayer-hydrate (BLH) is 
now attracting mounting attention due to the recent dis- 
covery of its superconductivity 1 . One of the interest- 
ing characteristic feature of this superconducting phase 
is that it is derived through a soft-chemical process for 
the parent material (Nao.7Co02). Namely, a part of the 
Na ions located between the C0O2 layers are extracted 
and H2O molecules are inserted for this region. Here, the 
role of H2 O is of great interest from both a chemical and 
a physical point of view, since it can be a key to under- 
stand the pairing mechanism, and furthermore, it may 
provide some guiding principles to synthesize other novel 
superconductors in layered transition metal oxides. 

It has been proposed from the beginning 1 that the 
main effect of inserting H2O molecules is expanding the 
interlayer distance (d) of the C0O2 layers. Namely, while 
d =5.4 A for the parent material, that of superconduct- 
ing BLH is expanded to d =9.8 A. Electronic structure 
calculations within the density functional theory 2,3 have 
indeed shown that the electronic structure of BLH is 
highly two-dimensional, whereas the c-axis dispersion of 
the band structure as well as the bonding-antibonding 
band splittings due to the interlayer couplings are not 
negligible for unhydrated NazCoC^. In fact, the two- 
dimensional electronic structure is favorable for uncon- 
ventional superconductivity where the pairing interaction 
has a characteristic structure in the momentum space 4,5 . 

On the other hand, to provide a deep insight into the 
role of the H2O molecules, Takada et al. synthesized 
monolayer- hydrate (MLH) by partial extraction of H2O 
from BLH, and studied the magnetic properties 6,7 . One 
may expect that MLH is also superconducting because 
d =6.9 A, which is larger than that of La2Cu04 (=6.6 
A). However, they found that the superconductivity is 
completely suppressed, while the carrier density in the 
C0O2 layer is expected to be unchanged because the mag- 
netic susceptibilities of MLH and BLH have similar values 
above 20K. Thus Takada et al. concluded that the H2O 
molecules in BLH may play an important role for the 
superconductivity besides the separation of the CO2 lay- 



ers. In fact, in contrast to MLH where the H2O molecules 
and the Na ions are placed on the same plane, in BLH, 
the H 2 molecules are placed between the C0O2 layers 
and the Na ions. Therefore, they proposed that the H2O 
molecules in BLH may shield the (random) Coulomb po- 
tential of the Na ions. 

While intensive first-principles calculations for BLH 
have been done up to present 8-12 , those for MLH have yet 
to be carried out. The purpose of the present study is to 
investigate why BLH is superconducting but MLH is not 
by comparing the electronic structure of BLH and MLH, 
and to obtain some hints to understand the mechanism 
of superconductivity in BLH. We found that BLH and 
MLH have similar band structures, but a notable differ- 
ence can be seen in the a\ g band around the Fermi level 
in that while its dispersion in the z direction is negligibly 
small for BLH, it is of the order of 0.1 eV for MLH. 

Although various theoretical studies have been per- 
formed for the superconductivity in BLH, the mecha- 
nism is still controversial. Especially, clarifying which 
bands among the ti g bands around the Fermi level dom- 
inate the superconductivity is of great interest. A va- 
riety of scenarios have been proposed. In fact, while 
there have been several works which claim that the a\ g 
bands are important 13-19 or both a\ g and e' g should be 
considered 20 ' 21 , the present author and his collaborators 
have proposed that the pocket-like Fermi surfaces of the 
e' g bands are essential 22 . Here, the advantage of first prin- 
ciples calculations is that it can give some hints on which 
band is important and should be taken into account for 
the model calculation. The present result that the main 
difference between BLH and MLH is the dispersion of 
the ai g bands suggests that the a\ g band should not be 
ignored when we study superconductivity in BLH, and 
the three dimensional feature of the a\ g band may be the 
origin for the absence of superconductivity in MLH. 

In the present study, we perform a first-principles 
band structure calculation within the framework of the 
generalized gradient approximation (GGA) based on 
the density functional theory. In the calculations, the 
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exchange-correlation functional introduced by Perdew, 
Burke, and Wang 23 is adopted with the ultra-soft 
pseudopotential 24,25 in a separable form. The wave func- 
tions are expanded by plane waves up to a cut-off energy 
of 30.25 Ry, and k point meshes of 8 x 8 x 2 are used. The 
atomic configurations are optimized to minimize the total 
energy with the conjugate gradient scheme 26 . We take a 
\fia x \/3a x c supercell, in which six cobalt atoms are 
included. During the optimization, the lattice constant 
a and c(= 2d) are fixed to the experimental values, i.e., 
a=2.8169 (2.8345) A and c=19.645 (13.842) A for BLH 
(MLH) 6,27 . While a variety of charge and spin orderings 
are observed in Na^CoC^ for x >0.5, such orders have 
not been observed in BLH and MLH (x ~ 0.3). There- 
fore, we performed a density functional calculation for 
paramagnetic states without assuming any charge order- 
ings for the initial state of the self-consistent calculation. 

In this paper, we focus on the electronic structure 
forNa 1/3 Co02(H 2 0) 2/ 3 (BLH) and Na 1/3 Co0 2 (H 2 0) 1/3 
(MLH). On the other hand, recently, it was reported 
that there are not only Na ions but also oxonium ions 
((H 3 0) + ) between the Co0 2 layers, so that the formal 
valence Co should be +3.4, which is lower than +3.65 
as estimated from the Na content 28 . While this means 
that the Fermi level locates slightly higher than that of 
the present calculation, for simplicity, we neglect the ox- 
onium ions. However, as is discussed below, we expect 
that the shift of the Fermi level does not change the qual- 
itative tendency of the present result. 



and Nai/ 3 Co0 2 (H 2 0)i/ 3 , we first introduced a reference 
system where the Na ions and the H 2 O molecules is repre- 
sented by a corresponding electron doping into the Co0 2 
layers. Namely, we calculated the electronic structure of 
the Co0 2 layers with an increased number of electrons 
and a uniform positive background charge to make the 
system charge-neutral. This reference system represents 
a situation where the H 2 molecules completely shield 
the potentials of the Na ions. 

In Fig. 1(a) and (b), we show the band structure of 
the reference system where the lattice constants are set 
to those of BLH 27 and MLH 6 , respectively. Here, the 
eighteen bands comprising the i 2s manifold are split into 
six a\ g and twelve e' g bands. Note that only nine bands 
would be distinguishable if the two layers in the unit cell 
did not interact. We can see in Fig. 1 that while the 
band dispersion in (a) and (b) are almost identical to 
each other, for (b) a small band splitting can be seen in 
the a\ g bands just above the Fermi level. 

In Fig. 2, we plot p(z) — J |0k(r)| 2 G?xG?y, where <^>k 
is the Bloch wave function of the fc-point indicated with 
arrows in Fig. 1 on the a\ g and e' g . Here, only the lower 
half of the unit cell is shown, and the Co planes reside 
around z = and 0.5c. We can see that the main ampli- 
tude of p(z) is localized within the Co layers. For BLH, 
the peaks on the Co planes are well separated for both 
dig and e' g . For MLH, on the other hand, the a\ g state 
has finite amplitude there whereas the e' g state vanishes 
z ~ 0.25c. 
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FIG. 1. The band structure of C0O2, for which the effect 
of the Na and the H2O molecules is represented by increasing 
the number of electrons in the system. The energies are mea- 
sured from the Fermi level. The lattice constants are set to 
be those of (a) BLH and (b) MLH. A small band splitting can 
be seen for the ai g bands just above the Fermi level. Solid 
arrows point onto the ai 9 and e' g band at the fc-points for 
which the Bloch wave functions are shown in Fig. 2. 



Before the actual calculation for Nai/ 3 Co0 2 (H 2 0) 2 / 3 
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FIG. 2. Squared absolute value of the Bloch wave function 
at the fc-points indicated in Fig. 1 for the bonding ai g (solid 
line) and e' g band (dotted line). Only lower half of the unit 
cell is shown and the Co planes reside around z — and 0.5c. 
The inset is an enlargement of p < 0.25. 



These results suggest that the assumption that the 
effect of the Na ions and the H 2 molecules is irrele- 
vant can be invalid for the a\ g states in MLH. Especially, 
the effect of the random potential of Na can be serious 
for these states. In other words, it is suggested that (i) 
for unhydrated Naa;Co02 2 ' 3 , the c-axis dispersion in the 
band structure or the interlayer couplings are not negli- 
gible for both e' g and a\ g orbitals because d is not suffi- 
ciently large; (ii) for MLH, d is large enough to suppress 



the interlayer transfers or couplings for e' g but not for a± g 
orbitals, because the wave functions of the latter extends 
to the c-axis; (iii) for BLH, d is sufficiently large so that 
the whole t2 g manifold is highly two-dimensional. It is 
expected that this qualitative tendency would not change 
even if we increased the amount of electron doping 29 . 
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FIG. 3. Density of states for the band structures of Fig.l 



In Fig. 3, we plot the density of states for Fig. 1. We 
can see that (a) and (b) are similar to each other. Espe- 
cially, both have peaks just below the Fermi level, which 
are due to the van Hove singularity in the e' g bands 22 . On 
the other hand, the divergence at the upper band edge is 
weak for (b) due to the band splitting of the a\ g bands. 



the facing oxygen atoms of adjacent C0O2 layers, they 
locate the midpoints of the prism edge for MLH (Fig. 4 
in Ref 6 ). Since the atomic position of the H2O molecules 
are not given in Ref 6 , we started with several patterns 
of initial configurations, placing the oxygen atoms of 
the H2O molecules on the same plane as the Na atoms. 
While we show the result for the most stable structure 
(Fig. 4(b)), the following argument does not depend on 
the detail of the atomic configuration of H2O atoms 30 . 
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FIG. 5. A plot similar to Fig. 1 of the atomic configuration 
of Fig. 4. 




FIG. 4. The optimized atomic configuration of (a) 
Na 1/3 Co0 2 (H 2 0) 2/ 3 and (b) Na 1/3 Co0 2 (H20) 1/3 . 

Next, let us discuss the actual calculation for 
BLH and MLH, i.e., Na 1/3 Co02(H20)2/ 3 and 
Nai/3Co02(H 2 0)i/ 3 . For the former, starting with the 
atomic configuration by Jorgensen et al. (Fig. 8a of Ref. 
27), we performed a structure optimization. In Fig. 4(a), 
we show the side view of the optimized structure. As for 
the initial atomic configuration of the latter, we followed 
Takada et al. 6 As is stressed in Ref. 6 , it should be noted 
that the positions of the Na ions in MLH are different 
from those in BLH. Namely, while for BLH the Na atoms 
locate near the center of the trigonal prisms formed by 



In Fig. 5, we show the resulting band structure. The 
difference between Fig. 1 and Fig. 5 is whether the in- 
tercalated hydrate group is actually considered or not. 
By comparing Fig. I and Fig. 5, we can evaluate the 
effect of the potential of Na ions, the screening by the 
H2O molecules, etc. We can see that the effect for the 
e' g bands in BLH and MLH is similar to each other. Es- 
pecially, for MLH, the band splitting in the e' g bands 
just below the Fermi level is negligibly small, which con- 
trasts with the case of unhydrated Na x Co02, for which 
the band splitting is as large as 0.1 eV 2,3 . This result 
implies that inserting single layer of H2O molecules is 
sufficient to suppress the interlayer couplings for the e' g 
bands. 

On the other hand, the c-axis dispersion and the band 
splitting of the a\ g bands is of the order of 0.1 eV for 
MLH, while those for BLH are negligibly small. This 
suggests that inserting a single layer of H2O molecules 
is not sufficient to separate the a\ g states in each C0O2 
layer since they extend along the c-axis. Therefore, we 
expect for MLH that while the e' g bands are rather in- 
sensitive to the random potential of the Na ions, the a\ g 
bands feel them seriously. 

In Fig. 6, we plot p for the fc-points on the bonding 
dig and e' g band which are indicated in Fig. 5. The 
c-axis dispersion is predominantly determined by the di- 
rect overlap between the tails of the adjacent main peaks 



rather than the amplitude on the H20 molecules (minor 
middle peak). We can see that the overlap is finite for 
the dig band of MLH but negligibly small for the other 
cases. (Note that for the e' g band of MLH, p vanishes 
between the main peak and the small peak around the 
oxygen site in the H2O molecule.) 
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FIG. 6. Same as Fig. 2, but now for the atomic configura- 
tion of Fig. 4. 



should not be ignored when we study superconductivity 
in Na x Co0 2 ■ yH 2 0. 

Recently, a systematic angle-resolved photoemission 
study has been performed for wide range of Na 
concentrations 31 . They concluded that the Fermi sur- 
face consists of the a-y g band for x =0.3, which suggests 
that the a\ g band is important for the superconductivity. 
This is consistent with the present conclusion. 

The GGA calculation was performed with TAPP 
(Tokyo ab-initio program package), for which we received 
technical advices from Y. Suwa. We benefitted from 
fruitful discussions with K. Kuroki. We would like to 
thank M.A. Korotin for providing us a program of the 
tetrahedron method and K. Held for his critical reading 
of the manuscript. This research was supported by the 
Alexander von Humboldt Foundation. 



In Fig. 7, we show the density of states for the band 
structure shown in Fig. 5. While the overall qualitative 
feature of Fig. 7 (a) and (b) is similar to each other, the 
divergence at the upper band-edge for BLH is completely 
smeared out for MLH. 
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FIG. 7. A plot similar to Fig. 3, for the atomic configura- 
tion of Fig. 4. 

In conclusion, we found that BLH and MLH have quite 
similar band structures, but there is a notable difference 
in the dispersion of the a\ g band around the Fermi level. 
Namely, it is negligibly small in the z direction for the 
former, it is an order of 0.1 eV for the latter. This sug- 
gests that although a single layer of H2O molecules is 
sufficient to suppress the interlayer coupling for the e' g 
bands, it is not sufficient for the a\ g bands. 

If we assume that the three dimensional feature of the 
dig band can destroy superconductivity, we can under- 
stand why BLH is superconducting but MLH is not. In 
this sense, the present result suggests that the a\ g bands 
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